IF 



4^. 



IN THF, TTNTT F.n STATFS PA^^NT ANn TR AnFTV/TARKr OFFirF 

In re Patent Application of: 

SKLARe^«/, 

Serial No.: 09/501,643 

Filed: February 10, 2000 

For: Flow Cytometry for High 
Throughput Screening 



Art Unit: 1641 
Examiner: Gabel, G. 

Docket No: UNME-0070-1 



Director of the U.S. Patent and Trademark Office 
Washington, D.C. 20231 



REQUEST FOR PRIORITY 
TnVDF.R35U.S.C.§12n 



Sir: 

In the matter of the above-captioned application for a United States patent, notice is hereby 
given that the AppUcant claims the priority date of November 9, 1999, the filing date of the 
Provisional Patent Apphcation No. 60/1 56,946 under 35 U.S.C. § 120. 



Respectfully submitted. 




Jagtiani & Associates 

10379-B Democracy Lane 
Fairfax, Virginia 22030 
(703) 591-2664 

June 21, 2001 



06/65/2001 14:03 5(35273|B|5 UIM CVTOMETRV FC^GE 



A F8ow Injection FSow Cytometry System 
for On-Line Monitoring of Bioreactors 



Buf Zhao, Arvlnd Natarajan, Friedrich Srtenc 

Department of Chemicsf engineering and Materials Science, Biological 
Process Technology Institute. University of Minnesota. 240 Gartner Lab, 
1479 Gartner Avenue. St. Paul MN 55108; telephone: (612)'624'9776; fax: 
612-625-1700; e-mail: fried^cbs.umn.edu 

f^eceived 27 February 19$$: accepted 9 September 1998 



Abstract: For direct and on-line study of the physiologi- 
cal states of cell cultures, a robust flow injection system 
has baen designed and interfaced with flow cytometry 
(FJ-FCMl. The core of the flow injection system includes a 
microch amber designed for sample processing. The de- 
sign of this microchamber allows not only an accurate 
on-line dilution but also on-line cell fixation, staining, 
and washing. The flow injection part of the system was 
tested by monitoring the optical density of a growing 
£.co// culture on-line using a spectrophotometer The en- 
tire growth curve, from lag phase to stationary phase, 
was obtained with frequent sampling. The performance 
of the entire R-FCM system is demonstrated in three ap- 
plications. The first is the monitoring of green fluores- 
cent protein fluorophore formation kinetics in £co// by 
visualizing the fluorescence evolution after protein syn- 
thesis is inhibited. The data revealed a subpopulation of 
cells that do not become fluorescent. In addition, the data 
show that single*cell fluorescence is distributed over a 
wid range and that the fluorescent population contains 
cells that are capable of reaching significantly higher ex- 
pression levels than that indicated by the population av- 
erage. The second application is the detailed flow cyto- 
metric evaluation of the batch growth dynamics of E.coli 
expressing Gfp. The collected single-cell data visualize 
the batch growth phases and it is shown that a state of 
balanced growth is never reached by the culture. The 
third application is the determination of distribution of 
ONA content of a S. cerevfsiae population by automati- 
cally staining ceils using a DNA-spectfic stain. Reproduc- 
ibility of the on-line staining reaction shows that the sys- 
tem is not restricted to measuring the native properties 
of c lis; rather, a wider range of cellular components 
could be monitored after appropriate sample processing. 
The system is thus particularly useful because it operates 
automatically without direct operator supervision for ex* 
tended time periods. ® 1999 John WUey & Sons. Inc. Biotech- 
not Bloang 62: 609-617, 1999, 

Keywords: online flow cytometry; flow injection flow cy- 
tometry; automatic staining; bioreactor monitormg 



INTRODUCTION 

There is ample evidence that growing cells are heteroge- 
neous entities that differ from one atwiher in their physi- 
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ological states. Microbial heterogeneity raay arise due to 
phenotypic changes associated with ihe cell cycle (Dicn and 
Srienc, 1991; Sricnc and Dien, 1992), due to changes in the 
mtcroenvironraents of individual ceils (Dunlop and Ye, 
1990; Fooler and Dunlop. 1989), or due to mutations re- 
sulting in genotypic vanations in the population (Hall, 
1995). Thus, rates of growth-associated parameters, such as 
protein synthesis or substrate uptake, are distributed in the 
population. From a biochemical engineering perspective, 
these population dynamics have profound implications, be- 
cause the overall productivity of the microbial process de- 
pends on the contribution of each individual cell. One 
nnethod to estimate such dynamics is by the use of flow 
cjrtomctry. 

In a process environment, rapid, repeated, and long-term 
on-line manual analysis is usually impractical^ if not impos- 
sible. Therefore, a certain degree of automation is desirable, 
particularly for more complex analysis procedures such as 
flow cytometry. To date, a number of flow injection systems 
have been designed and widely used in microbial process 
control and automatic analysis (Munch ct al., 1992; Reed, 
1990; Ruzicka and Hansen, 1975, 1988), but most of them 
moTutor population-averaged properties. The concept of au- 
tomatic flow cytometric analysis was introduced by Omann 
et al. (1985), who constructed a sample introduction device 
fitted to a Becton Dickinson FACS cytometer. Kelley 
(1989) also developed a similar device. Pcnnings et al. 
(1987) developed a system based on continuous pumping of 
cells and reagents with a peristaltic pump. Although the 
capabilities of these designs are rather limited, they are 
pioneering designs in automated flow cytometry. Successful 
on-line flow cytometry has been demonstrated using a flow 
injection technique (Ruzicka and Lindberg, 1992). 

In this .study, we present a flow injection system inter- 
faced with a flow cytometer and a bioreactor to perform 
on-bne assessment of single-cell property distributions. The 
versatility and performance of this system are demonstrated 
in several preliminary examples that show the utility of the 
system as it provides detailed quantitative information on 
growing cell populations that cannot be obtained with any 
other existing method. 



© 1999 John Witev & Sons, Inc. 
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MATERIALS AMD WIETHOOS 

Cell Strain*. Growth M«d!um, Growth 
Conditions, and St«9ntng Conditions 

E. coll strain K12 wa$ used in the growth expenment in* 
volving on-line monitonng of optical density Cells were 
grown m 300 mL of 2XGYT complex medium (Sambrook 
et ai., 1989) in a l-L flask placed in a 30'C water bath and 
shaken at 225 rpm. The culture was aerated at I wm by a 
peristaltic pump through an airstonc (Cole- Parmer, Vernon 
Hilh, IL) immersed in the medium. coli BL21 ceils trans- 
formed with the plasmid pRSET/S65T (a kmd gift from Dr. 
R Y. Tsicn. Howard Hughes Medical Center. San Diego) 
were used in the expenmem to study the Gfp fluorophore 
formation kinetics and batch growth dynamics. The plasmid 
contains an arnpicillin-resistance marker gene» and the Gfp 
gene under the control of a T7 RNA polymerase promoter. 
Cells were grown in LB complex medium containing am- 
piciUin (lOO j*g/mL). Production of T7 RNA polymerase 
was regulated using an IPTG inducible lacZ promoter pre- 
sent in the host chromosome- IPTG (200 ^g/mL) was used 
lo induce Gfp expression, and chloramphenicol (30 ^$/mL) 
was used to inhibit the protein synthesis when needed (Sam- 
brook et al., 1989). 

Sacckaromyces cerevisiae strain YPH399a (AfATa. ad€2- 
101, l€u2AL ly52'S0, his3^200, trpl^63, ura3'52) cells 
were grown overnight in 3 mL YPD medium (bacto-yeasi 
extract 1% w/v, bacto-pepione 2% w/v, dextrose 2% w/v) at 
30''C and 225 rpm in a IS-mL polystyrene lest tube (Fal- 
con). Cells were diluted to a concentration of ca- I x 10^ 
cells in fresh medium^ and tlw tube was placed on ice. 
Samples were automatically withdrawn into the niicrocham- 
ber of the flow mjection system. Inside the microcharaber, 
samples were washed with ice-cold PBS (0.15 A/, pH 1,S\ 
fixed in ice-cold 70% cthanol, washed with ice-cold PBS, 
treated with chromatin deuaturation solution (O.IN HCl, 
0-5% w/v Triton X-lOO, 1.75% w/v NaCl), washed with 
ice-cold PBS, and stained with mithramycin A (Sigma, 30 
p.g/tnL in PBS, 2 mAf MgCl^)- Fixation and denaturation 
steps were 2 min in length, washes were I min m length, 
and cells were stained with mithramycin for 10 min. Fixa- 
tion, denaturation, and washing were performed by continu- 
ously pumping the appropriate reagent through the micro- 
chamber while mixing the suspended cells using a magnetic 
stir-bar. 



Fi-FCM System 

The equipment used to construct the flow injection system 
is listed in Table L The components were interfaced to a 
personal computer using DT2805 and DAS 160 1 data acqui- 
sition and system control boards through digital input/ 
output (Dl/O) and digi^^l-to-anaiog (O/A) subsystems. 
Labiech Notebook software (Laboratories Technology 
Corp., Wilmington, MA) was used to control both boards. 
Figure 1 shows a schematic overview of this FI-FCM sys- 



Table I. Lisi of cquipmem. 



Equipment 



Manttfac'urer 



L T€n-position 
switching valve. 
Model 

2. Two-way injection 
valve. Modal 

3. Three- way 
switching valve, 
Mo<tel 01367-72 

4. fcristallic pump, 
7520-50 

5. Syringe pump, 
pomp 22 

6. Peristaltic pump. 
7520-20, 7520-30 

7. Magnetic Stir 
Plate, McKlel IZO? 

8. LH fennentOT, 2L 



Valco InsinjmentJf 
Co. Inc., 
Houston. TX 

Vaico 



Cole- Partner, 
Vernon Hills. 

Cole- Pwmcr 



IL 



9. 



10. 



U 



12. 



[)T2805 data 
acquisition and 
5ys(em control 
bostfti 

DAS 1601 data 
acquisition aod 
system control 
Ijowrd 
Computer, 
Vectra/66 

Compuicr, 486/2 3c 



Harvard 

Instrumentation- 
Harvard. MA 

Cole-Parmer, IL 

Fi?hef Scientific, 
Piiisbttrgh, PA 

LH Fermenution, 
Haywurd, CA 

Data Trans latioo. 
Marlboro, MA 



Keithley Qata 
AcqQiSiti<m* 
Taunton, MA 

Hewlen Packard 

Gateway 2000, N, 
Sioux City, SD 



Control method "-^h: 

*- ' ■ " ' ■ ■ "j*i*V 

TTL output. Dl/O % 
port** 1.0X2805 

TTL output, Dl/O 
port #1, DT2905 

Relay driven. Dl/O 
port #0, DT2503 

Current output; D/A 
pon #K DT2805 
Board 

RS-232C 



Masterflcx speed 
control unit 



LH fermentor 
control unit 



Interfaced with 
flow cytometer 

Interfaced with FI 
lysietn 



tem, which consists of three subsystems: (i) sample deliv- 
ery; (ii) sample handling; and (iii) sample injection and 
analysis- 

A sample delivery loop transferred the cell culture from a 
biorcactor to the flow injection system. A static degassing 
unit was designed to release air bubbles trapped in the 
sample (Fig. 2a), and sample was continuously recirculated 
in this loop. During sampling penods, the cell culture with 
air bubbles was allowed to accumulate in the glass tube, and 
a weak vacuum was applied using a peristaltic pump. EHic to 
a combination of static hydraubc buoyancy force and the 
vacuum, air bubbles were rapidly eliminated from the me- 
dium through a 0.45* ^tm in-line filter. The degassed sample 
was then fed into the microchamber of the sample handling 
subsystem for further processing. The sample residence 
time in the delivery loop was mininuzed, because the envi- 
ronmental conditions (in particular, the aeration) in the mb- 
ing were not the same as diose in the bioreacior. However, 
care was taken to avoid shear-induced damage of the sample 
that might result from the use of very high flow rates in the 
tubing. Hence, Cole-Paimer Mastciflcx silicon tubing (size 
#13, i.d. 0.75 mm) and a Oow rate of 5 tnL/nun (0.19 m/s) 
were used, resulting in a residence time of 10.5 s. 

The sample handling consists of a ten-position switching 
valve (#1 , Table I) connected with a precise peristaltic pump 
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Figure I, Aa overview of the FI-FCM systern. The system con^bb of the 
foUowin^ parts : sample delivery ("S D/'): s^pje handling C"S.R"); 
jample injecuoo ("S.l."); sample analysis ("S.A *'); aud cojnputcr control 
("C.C/'). Sample from a bio«actor is first brought to a tnultipowiioti 
valve «uj degasaing is pcrfbnned. Up to 10 different sampler or reagents 
("REl/* etc.) can be 5elected and pumped into the system using the 
ten-po9itton valve connected with a piecise pcriftaJtic pump. A micro- 
chamber CMC") Connected to a two-way injection valve is uied to diliice 
sample^ or mix sample? i*ith different reagents. After trrotmeot, samples 
w injected into the flow cytometer in > pulse-free manner by the mobile 
phase (saline), which i$ driven by pressurized air. The flow injection sys- 
tem is controlled by a PC throMgh two date acquisition and syeiem control 
boards. Flow cytometric dau acquisitioii and analysea m perfbrtned by 
another PC. 

(#4, Table I) to 5ciecT up to ten different streams, and a 
two-way injection valve (#2. Table I) incorporated with a 
microcharnbcr to load and inject sainples (Fig- 2b»c), The 
key component in sample handling is the microchamber. 
which has been designed to allow on-line sample dilution 
and staining (Fig. 2d). The unit essentially represents a 
stirred-tank reactor with three ports that serve as inlets and 
outlets- Because dilution, staining processes, and other en- 
zymatic reactions are basically mixing processes, they can 
be easily carried out in the microchamber in a predictable 
manner. 

Ports A and B are directly connected to the microcham- 
ber, whereas port C is connected to the microchamber 
through an in-line filter. With this in-line filter, fluids can 
flow through the tnicroch^ber freely, but cells are rcuined 
inside. To load sample into the microchamber, outlet C is 
blocked (Fig. 2d), sample is pumped through microchamber 
from pon A to pon B. To perform on-line diludon of cell 
samples, water is pumped through the microchamber from 
port A to pon B at a predetermined flow rate, F, for a certain 
time, (, such rhar the sample is diluted by a factor D given 
by: 




Figure 2. (a) The ??atic degassing unit. Air bubbles trapped in the cell 
culture escape from die medium due to the buoyancy force enhanced by the 
weak vacuum provided by a pcri»talUc pump (PF2)- (b. c) Schematic of 
>ample handling jubsystcm. The center line of the ten-position switching 
valve i9 connected to a perijultic pump (PP). Up to ten different streams 
can be fed into the system using this pump. The syringe pump (SP) is 
specifically assigned to perform on-line dilution, and to flush the system 
between samples. A microchamber (MC) is connected m the two-way 
injection valve (d) Cross- scctionaJ and three-dimensional views of the 
microchamber. Two side connections (A and ft) allow fluids and cell 
psiTicles to flow through, and the vertical connection CO ia separated from 
the microchamber thrpugh a membrane that allows only fhiida to pass 
through, (c) Sample injection subsystem. 



(I) 

where V is the volume of the microchamber, and Co is the 
initial cell concentration. In practice, the volume term is 
modified to account for the dead volume of the connection 
tubing. To perform on-line staining, port C is opened, 
whereas pon B is blocked. Solutions such as PBS or ethanol 
are pumped through the microchamber from port A to pon 
C, so that cells inside the microchamber can be washed, 
fixed, and stained. Afiter samples are processed, they are 
injected into the flow cytomcter for artalysis. The flow cy- 
lomctcr used was an Ortho Cytofluorograf Us {Onto Diag- 
nostics Systems, Westwood, MA) equipped with an argon 
ion laser Hnnova 90-5, Coherent Inc.. Palo Alto, CA) and a 
Cicero data acquisition system (Cytomation, Fort ColUns, 
CO). The laser was operated at 485 nm (Gfp) or 457 nm 
(miihramycin). and 100-mW beam power. A band-pass fil- 
ter (525 ± 15 nm) and an OG530 long-pass filter were used 
to collect Gfp and raithramycin fluorescences, respectively. 
Dau were acquired in listmode, pulse area, linear, and loga- 
rithmic configurations. 

RESULTS 

System Performance 

To test the sample handling system, a growti) study of E.coli 
cells was carried out using a spectrophotometer as the de- 
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lection device instead of a flow cytoraeter. The purpose of 
the experiment was to te^t the reliability of the sampling, 
degassing, and the dilution subsyHtetns. A satisfactory per- 
formance of the system would yield the growth curve of a 
cell cuhurc with frequent sampling points over the entire 
growth period by automatically carrying out dilutions to 
keep absorbance readings in the linear range when higher 
cell densities were reached. 

The output transmittance signals of the spectrophotom- 
eter were acquired using the A/D subsystem of the DT2805 
board. A typical output signal is shown in Figure 3a. The 
signal reserables a gamma distribution due to the dispersion 
of the sample in the carrier stream. The satnplc concentra- 
tion was determined by integrating over the cnrire duration 
of the signal pulse. Because the dispersion of the cell culture 
within the carrier stream is not linearly related with the 
sample concentration, the recorded transnutiance signals are 
not directly proportional to the sample concentration. This 
relationship can be described by the Taylor dispersion (Tay- 
lor» 1953) as follows: 



<2) 



where C(f) is the concentration of the sample at the position 
of the detector at time r. Af is the total amount of the sample 
injected into the carrier stream* E is the dispersion coeffi- 
cient. R is the diameter of the tube, v is the velocity of the 
stream, and Z is the distance between the detector and the 
point of the injection. Transmittance [V{t)\ is related to cell 
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Figure 3. (o) Typical 5p«ctrQphptomcter output transrmttance signal The 
signal w«g acquired using the A/D subsystem of ihc DT2S05 data acqui- 
5Uton board. (b> The range ovrr which Beer-Lambert'? law 15 valid wns 
determined by diluting a dense culture, ami nte^ring ab&orbaiice at 600 
nm at each diltttiort. (c) £. coU growth curve. £. coli BJ-21 cclJi were 
grown in 2XGYT medium ai 30'C, and the optical density of tb« cell 
culture waj mqititoi^d on-line every iO min for 22 h. The biomas? doubling 
tinie was estimated as 1.77 h. (d) On-line dilution. Qn-ttne dilution started 
automatically 1 1 h after the st^rt of the experiment, and ihc dilution factor 
increased to 27 over the course of the cxpeTimcnt- 



conccntration [C{/)1 according to Beet-Lambert '.s law 
CZubay, 1993): 



(3) 



where K is molar absorption coefficient, and d is the path 
length. Combining Eqs. (2) and (3). the following expres- 
sion can be derived to relate the amount of sample to the 
measured transmittance signal: 




r 

Jo 



(1 - \ogV)dr ^ M 



(4) 



Eq. (4) is vahd over the same linear range described by 
Beer-LambcTt 3 law. This was determined to be the case for 
sample absorbances of <2.0 (Fig. 3b), Therefore, on-line 
dilutions were necessary for dense samples, A zero-order 
algorithm was used to calculate the dilution factor for the 
next sample (Z?, +1) based on the current concentration mea- 
surement (C() and the it>itial sample concentration (Cq). The 
algorithm: 

A>i=C/Co (5) 

is valid when the frequency of sampling is greater rjian the 
frequency of cell division. The lag phase, exponential 
growth phase, and the stationary phase can be clearly rec- 
ognized in the growth curve thus obtained (Figure 3c), On- 
line measurements showed very few variations associated - 
with measurement errors, and measurements were thus 
stable and consistent. The on-line dilution started automati- 
cally 1 1 h after the start of the growth expetiraent, axtd 
dilution factors increased from I to 27 over the span of the 
experiment (Figure 3d)- The degassing unit was able to 
efficiently remove air bubbles from the cell culture samplea 7 
even during the late stages of the experiment when the ceO 
culmre was very dense (OD --40). The two abnormal point* 
shown in Figure 3 a (at ca, 4 and 10 h) are probably due to 
air bubbles that likely originated from open pons on ih^:---^. 
ten-position switching vaive. This potential problem can be 7^; 
avoided by sealing unused ports. One can note that the r 
dilution factors calculated after the overshtxits in the mea^ ,V 
siirements did not diverge away, indicating that the zenv 
order algorithm is robust for this application. 

Reproducibility of On-tine Flow 
Cytometry Measurementd 

A Stable sample stream is very important to obtain accurate 
measurements on a flow cytometer. Hence, pressurized aif 
was used to drive sample from the tnicTOchamber to 
flow cytometer (Fig. 2e). Most FI systems, including sonfti 
Fl-FCM systems, use a mechanical synnge pump (Blankest, 
stein et ai,, 1996; Lindbcrg et aK, 1993) Using our s) 
pump» it was found that the coefficient of variation of 
scatter distributions of calibration beads obtained wat 
higher than 10% with strong concomitant background . 
In contrast, the use of pressurized air to drive samples 
reproducible, noise-free light-scatter distributions of 
bration beads. 
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To test whether the OT^-Une measureTnents indeed yield 
the correct distributions, replicates of light-scattering inten- 
sity distributions of 5. cerevisiae ceils measured on-line and 
off-line were compared (Figure 4a). It can be seen that the 
distributions obtained from on-Unc measurements can be 
superimposed over those from off-line measurements. The 
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Figure 4. (a) Reliability. The lighr-scaHjer distributions cf on-line and 
ofMiTie roeasurcmcms are compared lo test the reliability of the on-Iin« 
mtajujrcments (b) Reproducibifity. The coefficient ^ variition of forw«rd- 
and right-*flei« Hglit -jcanef (FALS and RAl-5) distribucion* of unifonn 
calibrmion beodv were measured aver an extended time period to test (he 
reproducibility of the oo-Une me«vreineDte. (c> Robustness. The time 
profiles of event frequencie& of samples of votiious concentration? (400* 
700. and 1700 events per i«cood. dufhed ttne$> were compared with 4D 
ideal profile thai would be obtained in the abiience of sample dispersion 
(v?Ud lines). The "sampling period" (time window when samples are 
acquired) and ' 'flushing period" of the sajnpie uialy^i? cycle am tndicaTsed. 



Kolmogorov-SinirnoY (KS) statistical test (Neter et aL, 
1988) was pcriformed to verify that the distributions can be 
considered statistically identical with a confidence level of 
>0.99. Thus, on-line measiurcments give the same informa- 
tion as off-line measurements. 

Reproducibility of on-line measurements as a function of 
time was tested by comparing light-scattering intensity dis- 
tributions of uniform calibration beads (2 013 ± 0.025 \xn\, 
Duke Scientific) measured every 10 min over a span of 4 h. 
Variance in light-scattering distributions of these beads is a 
usefiil indicator of the optical alignment of the flow cytom- 
eicr. Heture, this measurement was used to estimate the 
stability of the on-!jne measuremems. It was observed that 
the coefficient of variation (CV) of the distributions fluctu- 
ated very little over 4 h of periodic sampling (Figure 4b). A 
statistical analysis confirmed that the CVs were identical 
with 95% confidence (/-test; Neter et al.^ 1988), confirming 
that the measuremems and the instrument settings were not 
subject to a drift over lime. 

The flow rate of samples that can be applied in a flow 
cytomcter depends upon the sample concentration. The ar- 
rival of ceils of a given concentration at the analysis point in 
the flow cytoraeter is governed by Poisson statistics. Hence, 
sample concentrations have to be adjusted to ensure that the 
probability of two or more cells being analyzed simuUa- 
nequsly is minimized. Here, a threshold of 1 % was set for 
coincident events; that is, samples whose probability of co- 
incidence of two or more cells exceeded i% of (he prob- 
ability of single-event occurrence were diluted prior to 
analysis. The lower end of sample concentrations waw de- 
termined by the minimum number of events desired for 
acquisition and the time available for sample collection, in 
our system, the working range of sample concentrations 
resulted in event flow rates between 100 and 1000 events 
per second. To examine whether samples whose concentra- 
tions varied in this working range were processed similarly, 
their time profiles (i.e., number of events as a function of 
time) were analyzed after sample injection for a range of 
concentrations. Under ideal circumstances with no sample 
dispersion, these profiles would be described by exponen- 
tially decaying functions with the specific rate of decrease 
equal to the dilution rate of samples out of the microcham- 
bcT. The time profiles obtained experimentally were nor- 
malised to unit area to yield event-frequency profiles for 
ease of compari.son, and they were found to be identical 
(Fig. 4c). This observation has several important implica- 
tions- First, the window between the time of injection of 
sample and the dmc of analysis of sample was independent 
of the actual sample concentration. Second, the relative cell 
ntmiber density of the sample can be detennined as the 
product of the number of events counted within a fixed time 
window (e.g., the sampling period'*; Fig. 4c) and the fac- 
tor of dilution. Thjxd. the actual proj5les obtained wei^ com- 
pared with the itieal profiles just described, and they were 
found to be very similar. After reaching a maximum less 
than 10 s after the sample front reaches the analysis point in 
the flow cylometer, event frequencies decreased^ as de- 
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$cnbed by the ideal profile. This implies that the micro- 
chamber can indeed be approximated an ideal stirrcd- 
laiik reactor. Finally, any residual sample is completely 
flushed away during the Hushing penod, and hence there is 
no cToss-contamination between successive samples. Tims, 
the system was robust in handling a range in ceil concen- 
trations, and yielded very reliable data. 

Monitorfng Cell Population Dynamics 

To simplify the process, wc initially measured native prop- 
erties of cells. These arc usually restricted to light-scatter 
intensity measurements. Whereas these parameters arc very 
useful in discriminating cell populations and, in some cases, 
in e.^iimatmg cell viability, they usually cannot be related 
directly to intrinsic quantitative biological properties of 
cells. Hence, we expressed the green fluorescent protein 
(Gfp) to obtain physiological data from which biological 
iiiformation can be directly inferred. Gfp is a naturally fluo- 
rescent protein that has been expressed in many heterolo- 
gous hosts (Cubitt et al., 1995), and is a very useful reponer 
for quantitative, noninvasive detection of single- cell gene 
expression (Natarajan et al,, 1998; Subramanian and Srienc, 
1996). Gfp expression was induced by the addition of IPTG 
to exponentially growing coU cells. The objective of the 
experiment was to examine the heterogeneity in levels of 
protein expression in single cells in an exponentially grow- 
ing population. Hence, protein synthesis was inhibited 5 
min after induction by addition of chloramphenicol. Thus, 
expression of Gfp in all cells was restricted to a **pulse'* of 
finite time. The time of induction is designated as lime "0" 
in the graph shown in Figure 5. The culture was sampled 
and analyzed every 10 min using online flow cytometry 
after induction and inhibition of protein synthesis. On-line 
dilution was not performed because inhibition of protein 
synthesis resulted in growth arrest and the cell concentration 
remained constant. Figure 5a shows the lime-e volution of 
the cellular fluorescence distribution. The population was 
initially nonfluorescent. After induction, the cells gradually 
increased in their fluorescence. The increase in the mean 
fluorescence value of the population (comparable to a popu- 
lation-averaged assay) is shown in Figure 5c (*'100%**). 
The kinetics of increase in mean fluorescence are sigmoidaK 
possible because nascent Gfp molecules have to fold, cy- 
clize, and oxidize the fluorophore scqueniialJy prior to turn- 
ing fluorescent. The increase in fluorescence saturates ap- 
proximately 3 h postinduction. However, a significant frac- 
tion of the population remains nonfluorescent (20%, Fig, 
5d). This subpopulation can be better visualized in the cy- 
togram of FALS vs. fluorescence (Fig. 5b). In addition, 
there is significant heterogeneity in expression among fluo- 
rescent cells, indicating variability in the capacity to pro- 
duce heterologous protein in the population. From a bio- 
tec hnoiogical process perspective, uniform maximal expres- 
sion of protein is desirable among all cells in the population. 
One should note that this maximum expression level cannot 
be inferred from the mean of the entire heterogeneous popu- 




PTgHW 5. (a) The evoluiion of nuorescertt populations. Orecn Huores- 
cencc intensity frequency distributions of the culture are plotted as a func- 
tion of tiOM. ih) The fluorescent population heterogeneity. A distinct popu* 
Ution of nonfluorcjcem cells *a3 observed, (c) The specific rate of in- 
crea.se of nuoresccnce of the top 5% and 10^ fluorc.iccnt fractions is 
greater than tfte mean specific rate of the entire population, (d) The frocuon 
of nonfluoreseeni cells was determined by comparing fluorescence inten- 
sity diMributions of s»mpic5 with those obtained from E. coli 0L21 cell* 
Qot expressing Gfp (ne^tive control). The firacdon decreased to a final 
value of ca- 2im. Tliw population probably represents plMmid-frec cellj. 



lation due to the presence of nonproducers and low-level 
producers. Rather, a single^cell assay of the kind described 
here is necessary. In the experiment just discussed, a choice 
of the top 10% or 5% of fluorescent cells indicates that this 
maximum capacity is three- or four-fold greater than the 
population average value (Fig. 5c). 

Monitoring Batch Growth Dynamics of £. coli 

To test the system over longer time periods, the FI-FCM 
system was used to monitor the growth dynamics of E. coli 
cells from lag phase to stationary phase. On-line dilution up 
to 700-fold was performed over the course of growth. E. 
coil BL21 cells transformed with plasmid pRSET/S65T 
were grown at 37^C in complex mediiun containing ampi- 
ciilin as a selection marker. Even though the culture was 
gro\m under noninducing conditions, uninduced "leaky" 
expression of Gfp was observed in the population. The cell 
culture was sampled and analyzed every 10 min, and for- 
ward-angle light scattering (FALS), nght-angle light scat- 
tering (RALS), and fluorescence intensity signals were col- 
lected. 

After a short lag phase of ca. 20 min. cells grew expo- 
nentially for 4 h. then entered stationary phase (Fig. 6a). The 
time evolution of the size (approximated by FALS intensity) 
distribution of growing E, coli culture from lag phase to 
stationary phase is shown (Fig. 6b). If the growth of cells 
was balanced, then frequency diatributions of various prop- 
erties such as cell si^e or protein content would be expected 
to be lime-invariant. In the early exponential growth phase. 
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FfauT* (a) The growth curve The relative number density of the sample was me^ured by counting the number of ceils aaalyicd dunng (he sampling 
p«iad and by 3cm«g the ttumto density of the fir,i sample to I. (b) The size distribution of growmg E^coH culmrc. The fonvard-angle l^»h^sc*«.r mten^.ty 
colkcted every 10 min to characterize the size distribution, (c) Meaii iDtcn«hrs of Ught scatter (corr^Utrd w,tb c«U size: fiUed tnangles) and 
fluorcacence (open diamonds) of the population changed over the cour« of cell growth. Dunng exponcnti^ growOi. cell size mcreascd while prweiD content 
decreased Tlie reverse trend w« ob5crvcd m the stationary phase, (d) The time e^olutiotl of the fluot^scence intenjjity distributions is repre^^cnted m 
isometric and contour plot*. The majority of the cells were nonfluorcsccnt during e^^pooemial growth. Celt* became incrra.mgly more fluorescent dunng 
late exponenttaJ/stalionary phase, untij, at the end of the ea^periment. two populations couW be disungusshed. 



the peak of the FALS distribution shifted to right (indicating 
ao increase in mean cell size), reached a maximum value, 
shifted back to lower mean intensity during late exponential 
growth. A similar trend is observed in the mean FALS 
signal intensity (Fig. 6c). Thus, growth was exponential, but 
not balanced- FALS intensity distributions were unchanged 
duriLRg the statiooaiy phase. The fluorescence intensity dis- 
tribution of the growing culture wa^ also monitored. Al- 
though Gfp synthesis was not induced, there was sufficient 
"leaky** expression of the protein for the fluorescence to be 
detectable, Intercsiingiy, heterogeneity of expression wxs 
observed even among uninduccd cells (Fig. 6d). Levels of 
'Meaky'' uninduced expression were seen to be inversely cor- 
related to the growth phase, as lag- and stationary-phase cells 
were more fluorescent than exponentially growing cells (Fig. 
6c)- Additionally, protein content distributions in the popu- 
lation also changed with time during the exponential growth 
period, funher indicating that growth was unba^lanced, a^ a 
time-invariant physiological state is never reached. 

Monitoring Cell Cycle Distribution of S, c^wisiM 

To demonstrate the staining capability of tlie H-FCM sys- 
tem, DNA contents of 5- t:ere\nsiae cells were measured 



on-line as described earlier- Even the application of a 
simple, single stain requires multiple preparatory steps dur- 
ing which cell samples might potentially be lost. To elimi- 
nate this possibility, a 0.2- pum in-line filter was used lo 
retain cells while they were treated with various chemicals. 
The flow through port C was analyzed for presence of cells 
after each step to ensure the integrity of the membrane. 
Membrane clogging, frequently a problem m many mem- 
brane systems, was virtually eliminated in our design be- 
cause of coexisting tangential flow between ports A and B, 
continuous vigorous stirring in the microchamber, and di- 
lution of cells to a concentration of ca lO^/mL prior to 
initiation of staining. Vigorous stirring also ensured that 
cells did not clump together during the fixation step. 

To demonstrate staining capability, replicates of a sample 
of 5. cerevisiae cells in late exponential phase were washed, 
fixed, treated with chromatin denaturation solution, stained 
with mithramycin (MI), and analyzed on the FCM. An asyn- 
chronous population would be expected to have a bimodal 
DMA distribution, where the two modes correspond to cells 
in Gl and G2 cell cycle phases, and this was observed (Fig. 
7b). Because the culture was in late exponential/early sta- 
tionary phase, a majority of the cells were observed in the 
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Fi^puTv 7. (b) A bivariate hiscograni of PNA comcnt and forward- angle 
Ught-scancr imensiiy (FALS) of A Uic expoi^otial 5. cerevijiae culture. 
The two subpopulations observed cwr«5pond tp cell? in the G I and G2 cell 
cycle phases, (c) DNA content dtauibuUM of three replicate samples of S. 
cenvisiae cells (solid !io€. Y2 wis) are compared with autofluorescence 
measurements of a fixed but unstained population (dashed line. Yl axis). 
DNA distributions ^ve^e obtained by staining with MI (see text). All dis- 
tributians are nonnalized to unit area for ease of comparison. AutoHuQ* 
resccnce distribution (shown on a slightly different scale) ia clearly uni- 
modal. and well separated from DNA distributions- (a) However, staining 
does iKJi affect light-sc^ittrr intensity nicasurrmcnis, and FaLS diatribn- 
tions of stained and unstained S.cerevisiae ceils are identical. 

Gl phase^ (Fig. 7c). DNA disiribPlions obtained were com- 
pared with autofluorescence measureincnts on fixed and de- 
natured cells that had not been stained. Autofluorescence 
distribution is clearly unimodal, and well separated from the 
stained distributions (Fig. 7c). In comparison, staining did 
not affect light-scatter intensity measurements^ hence FALS 
distributions of unstained and stained cells were identical 
(Fig. 7a). DNA distributions for three replicate samples 
were also found to be statistically identical, ensuring repro- 
ducibility of the staining process (KS test; Neter. 1988). 
Because the entire process is automated, it is easy to envi- 
sion the use of more complicated staining protocols to lag 
and quantiiate other cellular components also. 

DISCUSSION 

To determine the state of a growing cell culture, conven- 
tional monitoring systems typically evaluate measurable 
quantities of the abiotic phase from which the state of the 
culture is inferred Direct measurement of the physiological 
state is usually not possible because quantitative methods to 
rapidly assay the composition of the biomass are limited. 
Furthennore, the properties of individual cells arc different 
and distributed over a range of physiological states. There- 
fore, population average data provide only lijnited informa- 
tion, and accurate information must be sought at the single - 
cell level. 
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Fredrickson et al. (1967) and Ramkrishna ct al. (1968) 
proposed a general mathematical framework, known as 
population balance theory, to describe cell growth and its 
interaction with the environment (see also Ramkrishna, 
1979, 1985). In this theory, the state of mdividual cells is 
specified by the physiological state vector — a collection of 
state properties such as cell size, protein contenu or DNA 
content. The resulting population balance equations are es- 
.semially number balances on individual cells of a popula- 
tion, which keep track of, not only the generation and dis- 
appearance of cells, but also of the continuous change in the 
identity of cells due to physiological growth processes. This 
modeling framework (known as corpuscular and structured) 
most realistically represents the evolution of a heterogc* 
neous microbial population. Ver?!?on5 of this model have 
been used to extract growth parameters from single-cell 
property distributions obtained using flow cytometry (Dicn 
and Srienc, 1991; Kiomenaker and Srienc, 1991, I994a-c). 
However, these analyses were performed off-line for special 
cases. In general, it has been difficult to apply the popula- 
tion balance model to real microbial processes due to lack of 
experimental data that would permit identification of model 
parameters. Such data could possibly be generated by flow 
cytometry interfaced with appropriate instrumentation for 
on-line sample processing and analysis. 

Hence, we designed and consimcted a versatile and 
modular flow injection-flow cytometry (FI-FCM) system. 
To display the details of information that can be obtained, 
the system was used to monitor the heterogeneity of protein 
expression in an E.coU cell population. The results cleariy 
reveal the need for single-cell analysis, because cellular 
content of even a single component such as Gfp content was 
distributed in the population. Frequent, automated sampling 
permitted observation of smooth trends of cellular phenom- 
ena. These data could be used to generate accurate estimates 
of the kinetics of transient events. Furthermore, it has been 
shown that the system can be used to monitor long-term 
fermentation processes, during the course of which cell con- 
centrations would change significandy. Because the system 
does not require any operator intervention, it is particularly 
suitable for process monitoring and control applications. 
The data obtained can be used to determine not only specific 
rates of exponential growth, but also to detennine the de- 
tailed physiological state distribution of the cell population. 
Additionally, single-cell growth parameters can be evalu- 
ated. To better define the composition of single cells, the 
range of cellular components that can be monitored can be 
expanded by raking advantage of the design of our micro- 
chamber, which permits on-line cell fixation, and staining, 
as we have demonstrated. 

The examples presented demonstrate that the developed 
system represents a powerful tool for studying cell physi- 
ology at a level of detail that was not previously possible. 
Moreover, the Fl-FCM system will l>e advantageous for 
monitoring biorcactors in studies involving transient growth 
phenomena, estimation of cellular kinetics, or study of syn- 
chronous or cyclic trends that might not be observed with 
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infrequent sampling. The ^ysiem ihm has significant poten- 
tial in biotechnology for the monitoring of microtnal popu- 
lations, without the need for excessive manual intervention. 

This work was wpponed partially by vhc National Scteoce Foun- 
dation (BES-9708146) and by ^ ?ced gram from thr Biological 
J*roc«j Technology IriKituic, Univeraiiy of Minnesota. 
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